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f*"""- ' The well-known Giudice-Masiero mechanism explains the presence of a fx term 

of the order of the gravitino mass, but does not explain why the holomorphic mass 
term is absent in the superpotential. We discuss anomaly-free discrete symmetries 
which are both compatible with SU(5) unification of matter and the Giudice- 
^| ■ Masiero mechanism, i.e. forbid the fx term in the superpotential while allowing 

the necessary Kahler potential term. We find that these are Zj^ symmetries 
with the following properties: (i) M is a multiple of four; (ii) the Higgs bilinear 
^ ■ H u Hd transforms trivially; (iii) the superspace coordinate 6 has charge M/4 and, 

$_i . accordingly, the superpotential has charge M/2; (iv) dimension five proton decay 

operators are automatically absent. All Zjy- symmetries are anomaly-free due to 
a non-trivial transformation of a Green-Schwarz axion, and, as a consequence, 
a holomorphic fi term appears at the non-perturbative level. There is a unique 
symmetry that is consistent with the Weinberg operator while there is a class of 
Tiff symmetries which explain suppressed Dirac neutrino masses. 
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1 Motivation 



The minimal supersymmetric standard model (MSSM) is a very appealing extension of 
the standard model of particle physics. Supersymmetry promises to stabilize the elec- 
troweak scale against radiative corrections. The structure of matter hints at unification, 
and the attractive picture of precision gauge unification pQ enabled by supersymmetry 
introduces the scale of grand unification Mqut = a few x 10 16 GeV. The MSSM also 
provides a compelling dark matter candidate. 

On the other hand, the MSSM has various problems. Usually the MSSM comes with 
matter or R parity [213] which eliminates the most troublesome baryon number violating 
interactions, and ensures the stability of the aforementioned dark matter particle. Yet, 
even after imposing matter parity, there are certain serious shortcomings. One of them 
is the so-called u fi problem" which consists in the question why the holomorphic mass 
term for the Higgs bilinear is of the order of the electroweak scale. In addition, there is 
the dimension five proton decay problem [3H5] (cf. also [6]). 

It is hence clear that the MSSM requires additional ingredients beyond matter parity. 
In this study we analyze anomaly-free discrete symmetries which forbid the \i term. As 
we shall demonstrate, requiring that the symmetries be compatible with the Giudice- 
Masiero solution [7J to the /i problem and SU(5) leads to very restricted classes of 
solutions, depending on whether neutrinos are Majorana or Dirac particles. In the first 
case, the solution is unique and even compatible with SO (10) while in the second case 
the smallness of the Dirac neutrino Yukawa coupling can be related to the suppression 
of the /i parameter. 



2 Naturally suppressed /x term and Dirac neutrino 
Yukawa couplings from anomaly— free symmetries 

We start by reviewing the explanations of a suppressed fi term through Kahler potential 
terms in Section I2TT1 Next, we discuss anomaly constraints in Section I2T21 In Section [231 
we comment on proton decay operators and study settings with SO(10) relations in 
Section 12.41 Then, we discuss the appearance of a suppressed holomorphic fi term and 
Dirac Yukawa couplings in Section 12.51 and Section \2.6\ respectively, and give a short 
recap in Section I2TT1 



2.1 Giudice— Masiero mechanism 

The famous Giudice-Masiero mechanism [7J provides a solution to the /x problem in the 
MSSM. Giudice and Masiero pointed out that in supergravity an effective holomorphic 
H u Hd bilinear, i.e. an effective /i term, can arise from the (non-holomorphic) Kahler 
potential term 

K d k HuHd ^-H u H d + h.c. (2.1) 
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Here X is the (spurion) field that breaks supersymmetry and kn u H d arid Mp denote a 
coefficient and the Planck scale, respectively. Inserting the F term vacuum expectation 
value (VEV) F x of X leads to an effective superpotential term 

^cff ~ -77- H u H d =: fi eS H u H d , (2.2) 

with /j, cS of the order of the gravitino mass m 3 / 2 , which sets the size of soft superpartner 
masses in gravity mediation. 

However, for the Giudice-Masiero mechanism to work, the holomorphic superpoten- 
tial term \i H u Hd needs to be absent in the first place, or better forbidden by a symmetry. 
As it turns out, symmetries that can forbid the \i term are rather constrained. It has 
been shown [5] that, if one requires the symmetry to be anomaly-free and to commute 
with SU(5) (in the matter sector), it has to be an R symmetry (cf. the similar discussion 
in j9]). As shown by Chamseddine and Dreiner [10], in the MSSM gauged anomaly-free 
continuous R symmetries are not available. On the other hand, there are strong argu- 
ments against global symmetries (cf. [TT] for a recent discussion). We are hence led to 
the conclusion that the symmetry needs to be discrete. In what follows, we therefore 
will only consider anomaly-free discrete R symmetries. Specifically, we will look at one 
particular generator which forbids the fi term. This generator will generate an Abelian 
discrete R symmetry of order M, i.e. a 1^ symmetry. 

One can actually narrow down the potential symmetries even further. Suppose we 
seek to generate an effective \x term, Equation f)2.2p . from the Kahler potential term. 
Here we assume that X is the field that breaks supersymmetry and generates gaugino 
masses. Then its F component has to have minus the R charge of the superpotential. 
One way to see this is by recalling that gaugino masses get induced by the operator 
f d 2 9 X W a W a (with 9 and W a denoting the superspace coordinate and the gauge mul- 
tiplets, respectively). Since the superpotential R charge q^ equals twice the R charge 
of 9, qo, and the lowest components of W a (i.e. the gauginos) carry R charge q#, the 
X superfield needs to be inert under the (discrete) R symmetry. Therefore, the Higgs 
bilinear H u Hd needs to be neutral as well. Altogether we have found that an anomaly- 
free and SU(5) compatible symmetry that forbids the \x term in the MSSM has to be 
discrete, and under this symmetry, 

e 27Ti M9, (2.3a) 

e 27ri T^ where q^ = 2q g , (2.3b) 

X , (2.3c) 

H u H d . (2.3d) 

Here and throughout this study we normalize the discrete charges to be integer, i.e. 
qe e Z. 

It is immediately clear that such a symmetry allows effective superpotential terms of 
the form 

W D Cn ^l HuHd ' (2 ' 4) 



9 


->■ 


W 


— >■ 


X 


— >■ 


H u H d 


->■ 
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where Q (with R charge q%/) denotes the superpotential of some 'hidden sector'. As 
usual, a non-trivial VEV of Q is required to cancel the vacuum energy. This VEV 
will break the R symmetry, but the breaking is hierarchically small, i.e. of the order of 
the gravitino mass 7713/2 (cf. the discussion in [T2] ). That means that, apart from the 
Giudice-Masiero contribution, one would expect to have a holomorphic ('Kim-Nilles 
type' [13]) contribution to the \x parameter of the right size. 



2.2 Anomaly constraints 

Up to now we have only used the fact derived in [5] that SU(5)-compatible and anomaly- 
free non-i? symmetries cannot forbid the fi term. Now we discuss anomaly constraints 
on lijh symmetries. These constraints have been re-derived recently in [8]. However, 
there only the special case qg = 1 has been considered, which is too strong a require- 
ment. To see this, consider a symmetry, for which there are two different non-trivial 
possibilities for the superspace charge, qg = 1 and qg = 2. At first glance, one may think 
that one may rewrite the qg Tif x Ti2 symmetry. This is not the case since 

2 and 8 are not coprimeQ The generalization of the anomaly coefficients to arbitrary qg 
is straightforward and deferred to appendix |A] 

After summarizing the relevant anomaly coefficients for the MSSM in Section 12.2.11 
we explain in Section 12.2.21 why 'anomaly universality' must be imposed in models in 
which the SM gauge group is unified into a simple gauge group. We then proceed by 
verifying the consistency with anomaly matching in Section 12.2.31 and show that only 
R symmetries can forbid the \x term in Section 12.2.41 Finally, we derive constraints on 
the order M in Section 12.2.51 and comment on the (ir)relevance of the universality of the 
mixed hypercharge anomaly in Section 12.2.61 



2.2.1 anomaly coefficients in the MSSM 

In the case of the MSSM the anomaly coefficients := A 

^SU(2) L -SU(2) L -Z« and A l '■= Aj(l) y -U(l)y-Z« rea d 



SU(3)c-SU(3) c -S« 5 A 2 



A R 



A* 



At 



R 



1 

9=1 

1 3 1 

9=1 

|E( 3 "?o+«f 5 ) + ? ri 



5g 



9=1 



(qH u + qH d ) ~ Uqg 



(2.5a) 
(2.5b) 
(2.5c) 



Here, q{ and q^ denote the SU(5)-universal R charges of the MSSM superfields (Q 9 , U 9 , E 9 



and (D 9 ,L 9 ), respectively, and g represents the flavor index. Accordingly, matter 



1 A simple way of seeing this is to recall that all elements of TL^ x Z2 have the property that taking 
them to the fourth power yields identity, which is obviously not the case for 7Lf. 



3 



fermions and Higgsinos have charges q — qe while gauginos have charge q$. 



2.2.2 Anomaly universality and discrete Green— Schwarz mechanism 

If the standard model gauge group is to be unified into SU(5) or SO (10), a necessary- 



condition for anomaly cancellation is the universality (cf. the discussion in appendix IB. 3|) 
Af = A* = Af = p mod 77 . (2.6) 

Here we introduce 

f M/2, if M even, 
' ' \ M , if M odd . 1 ' 

p is a constant which indicates whether or not a Green-Schwarz (GS) mechanism [2] is 
at work. Specifically, p is related to the discrete shift of the GS axion (see Equation ( 1B.9|) 



in appendix IB.2j) . p = means that the symmetry is anomaly-free in the conventional 
sense, i.e. without GS mechanism. 

At this point, we would like to comment on certain important properties of the 
Green-Schwarz mechanism and its discrete version as there seems to be some confusion 
in the literature: 

1. Although the GS mechanism plays a prominent role in string theory, it does not 
rely on strings. In fact, as shown in appendix IB.2[ it can entirely be understood 
in (the path integral formulation of) quantum field theory. 

2. Unlike in the continuous case, for discrete symmetries the transformation of the 
axion is only fixed modulo r\. It will be interesting to see whether this ambiguity 
can be fixed somehow, e.g. in explicit string-derived models. 

3. In the continuous case, the axion has to be massless for the shift symmetry to 
be a symmetry of the Lagrangean. That is, the axion potential needs to be flat. 
By contrast, in the discrete case the potential is only required to be periodic, 
i.e. invariant under the discrete shift, Equation (IB. 9}) . Therefore the axion may 
have a non-trivial mass prior to the breakdown of the symmetry. This is, in 
a way, somewhat surprising as it means that a massive (and bosonic) state can 
contribute to an anomaly. Of course, in both cases the symmetry will be broken 
(spontaneously) once the axion a acquires its VEV. 

2.2.3 Anomaly matching 



It is instructive to use 't Hooft anomaly matching [T5] (see [H] for discrete anomaly 

-com 

^SU(5) 2 



matching) in order to constrain the properties of anomaly-free GUT-compatible Wif 



symmetries. At the SU(5) level, there is only one anomaly coefficient ^4su(5) 2 -z^ f > which 
we can split into three parts, 
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The first term contains the contribution of matter and is given by 



^u(5)U« = \ E ( 3 ^o + 4) - 6?* • (2-9) 

9=1 



Here, we used Equation (lA.4aj) with Dynkin indices £(5) = §-£(10) = \. The second 
term in (12.81) . ^gu(5) 2 -z fl > denotes the contributions of additional fields, e.g. the SM and 
SU(5) breaking Higgs. Finally, the last term in ( 12 .8p represents the gaugino contribution 
for SU(5). Yet, by considering the SU(3)c and SU(2)l subgroups of SU(5), one can 
introduce two anomaly coefficients A^^L wr — ^cn^ w « at the GUT level, 

^su(3)2,-z« = ^su(3)2-z« + ^su(3)2-s« + 3 * + 2 ' 2 ' 2 ' qe ' (2.10a) 
^su(2)2-s« = ^su(2)£-s« + ^su(2)2-s« + + 2 ' 2 ' 3 ' qe ' (2.10b) 

where we artificially split the gaugino contributions into those from the adjoint represen- 
tations of SU(2) L or SU(3)c, respectively, and in those coming from the extra gauginos 
in the (3, 2)_ 5 / 6 © (3, 2) 5 / 6 representation. Assume now there is some (unspecified) mech- 
anism that breaks the GUT symmetry down to the SM symmetry, and thus removes the 
extra gauginos, while leaving 2^ unbrokenJl Then, the coefficients 

4 SU(5) broken _ ,«SU(5) _ 9 (OA'i^ 

^su(3)2-s« ~ ^svwi-z* zqe ' l^.liaj 

4 SU(5) broken _ ^SU(5) _ r, fOIlM 

/1 SU(2)2- Z fl - /i SU(2)2-Z« ^« l^.llDj 

cannot be equal, i.e. the anomaly coefficients cannot be universal, unless there are split 
multiplets contributing to A^™ )2 _ a H ( wnere we use ^su(3)2 = ^su(2)2-s« )• That 
is, 't Hooft anomaly matching for (discrete) R symmetries implies the presence of split 
multiplets below the GUT scale. 



2.2.4 Only R symmetries can forbid the /i. term 

Given that SM matter furnishes complete SU(5) representations and the attractive pic- 
ture of MSSM gauge unification, arguably the most plausible candidates for such split 
multiplets are the Higgs fields. Requiring that the Higgs fields cancel the mismatch of 
gaugino contributions to the anomalies, we obtain 

2 (Qh u + QH d - 2<? e ) = qe mod rj , (2.12) 

2 If one is to obtain the exact MSSM spectrum after GUT breaking, this mechanism cannot be 
spontaneous symmetry breaking in four dimensions [17] . On the other hand, extra dimensions, especially 
in the framework of heterotic orbifolds, naturally can give discrete R symmetries as remnants of higher 
dimensional Lorentz symmetry, see e.g. |18j . 
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implying 



<lH u + qH d = ^qe mod 2?7 (2.13a) 
= 2q w mod 2rj (2.13b) 
^ q w mod M for q w ^ mod M . (2.13c) 

Therefore, to.oto.-R symmetries with qg = = cannot forbid the fi term. But in case 
of non-trivial symmetries (i.e. M > 3) the term will always be forbidden, as it 
should be, since only chiral contributions can 'repair' the gaugino mismatch. 

A remark is in order to show that Zf 7 with q^ = mod M is not an R symmetry. 
In the case qw = mod M we find two solutions for qe\ either qg = 0, such that the 
symmetry is clearly non-i?, or (for M even) qg = M/2. However, since the transfor- 
mation 9 i — y —9 and ^ h- > — for all fermions \1/ is always a symmetry, one can shift 
the charges by M/2 such that again q e = 0. Hence, 1^ with = mod M is 
equivalent to a non-i? symmetry |19j . 



2.2.5 Constraints on the order M 

Using Equation (12.131) and assuming a Giudice-Masiero-like mechanism such that qn u + 
qu d = mod M from (12.3dj) . we obtain 

2q w = mod M , (2.14) 

which implies, given the freedom to choose q^ between and M — 1, that the only non- 
trivial solution for even M is qy/ = M/2. For odd M there is no non-trivial solution. 
Since the superpotential charge is given by qw = 2qg, the order M has to be divisible 
by 4. Hence we can focus on 

M = 4 x integer and q e = M/4 (2.15) 

in the rest of our discussion. 



2.2.6 No additional condition from 



Subtracting A% from yields 



(qH u + qH d ) - llqe 



+ 3qg = mod rj 



(2.16) 



Since M is even and qn u + qu d = mod M by (12. 3d}) , this equation is equivalent to 
3kM + 2(15 -33) q e = UM (2.17) 

with some integers k and I. That is, 

36 qg = [3k-5£}M = 7L ■ M . (2.18) 

For a given order M, this relation constrains qg. However, we know already from our 
discussion below Equation ( I2.14p that M needs to be an integral multiple of 4, such that 
( I2.18P does not lead to an additional constraint. 
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2.3 Family— independent symmetries and proton decay 

Let us now assume further that the discrete symmetry be Abelian, i.e. of TL^ type, with 
family-independent charges. Assuming the presence of Yukawa couplings, symmetries 
with the above properties have automatically the virtue of solving the dimension five 
proton decay problem of the MSSM, as we will see in the following. 

The requirement that up- and down-type Yukawa couplings be allowed, 

2gio + = qw mod M , (2.19a) 
Qio + Q5 + QH d = qw mod M , (2.19b) 

implies 

3<?io + % + Qh u + qn d = ^qw mod M . (2.20) 
Imposing (12. 3d}) gives 



3<?io + <?5 = Iqw mod M ^ q w mod M , (2.21) 

for an R symmetry (i.e. for qw ^ mod M), showing that the troublesome dimension 
five operators 10 10 10 5 are automatically forbidden whenever the Yukawa couplings 
are allowed. In [8] the same conclusion was obtained from anomaly cancellation. 
Recalling further that 2qw = Aqg = M leads us to the conclusion that 

qg = - 3g 10 mod M . (2.22) 

This means that the contributions of matter fields to the anomaly coefficients ( 12. 5 p 
vanish, and that the universal anomaly coefficients are simply given by 

Af = p = q e mod M/2 (2.23) 

for 1 < i < 3. 

Next, one can also discuss proton decay originating from the dimension four operator 
10 5 5. This operator has R charge 

qw + 2<?5 = - 5gi mod M . (2.24) 

Hence, these operators are also forbidden if — 5gio ^ qw mod M, or equivalently 
lOgio ^ k M with k odd. 

2.4 Imposing SO(10) relations 

Let us now comment on the special case that the symmetry commutes with SO(10) 
for the matter fields, i.e. q 10 = q$ = qw- The requirements that the up-type and 
down-type quark Yukawa couplings be allowed imply that qn u = qjj d =: qu (modM). 
Furthermore, from the anomaly universality condition ( I2.13bj) we find qu = qw mod r\. 



In the following, we consider two cases: in case (i) we demand in addition the Weinberg 
neutrino mass operator, and in case (ii) a Giudice-Masiero-like mechanism. 
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(i) If we require the Weinberg neutrino mass operator, i.e. 2qie + 2qn = qw mod M, 
we find M = 4m, m G N and 



m , q-g/ = 2m , q H = and q 16 = m . 



(2.25) 



This symmetry automatically allows for the Giudice-Masiero term and the universal 
anomaly coefficients Af = m ^ indicate a discrete GS mechanism. The simplest 
case m = 1 is the Tif symmetry discussed in [2U1I2T]- All other cases are just trivial 
extensions as long as one considers the MSSM states only. Of course, if additional states 
are introduced, they can have 1if m charges in such a way that one cannot reduce it to 
Tif. Another version of the uniqueness proof of Wif can be found in [8j. However, the 
analysis in [5] assumed that qg — 1. Here we show that uniqueness also survives the 
generalization to general qg ^ 1. 

(ii) If we do not require the Weinberg neutrino mass operator but a Giudice-Masiero- 
like mechanism, i.e. 2qn = mod M, there are two cases: both cases have M = 4m, 
m G N, qg = m and q^/ = 2m. In addition, in the first case we get qu = as discussed 
above in case (i), and in the second one we find qu = M/2 = 2m and qie = 21 m with 
t G Tj. However, this choice forbids the Weinberg neutrino mass operator. 

2.5 Non— perturbative holomorphic \x term 

If the above discrete R symmetry appears anomalous, i.e. if anomaly freedom is due 
to a GS mechanism (see appendix [B] for a discussion of its discrete variant), then such 
holomorphic contributions will appear as arising at the non-perturbative level [8"|l2"Tj. 
To see this, recall that the superfield S containing the axion a, i.e. S\g=o = s + ia, needs 
to enter the gauge-kinetic function, or, in other words, Jzf D f d 2 9 fs S W a W a (with 
some coefficient fs). Non-invariant terms in the superpotential can be made invariant 
by multiplying them by e~ bS with appropriate b. As s controls 1/g 2 such terms go 
like e~ b ^ , i.e. have the form of instanton contributions. This then fits nicely into the 
scheme of dynamical supersymmetry breaking [22] (see also the more recent discussion 
on "retrofitting" [23]), where the scale for supersymmetry breaking is set by a gaugino 
condensate [24] . or a more complicated dynamical term (see e.g. [25] for a review of 
simple models). 

2.6 Small Dirac neutrino Yukawa couplings 

By relating them to supersymmetry breaking one may explain suppressed neutrino Dirac 
Yukawa couplings [2&H28] . That is, similarly to the \i term, one can get effective Dirac 
neutrino Yukawa couplings from the Kahler potential terms 



K D 




(2.26a) 



as well as 



K D 




v + h.c 



(2.26b) 



S 



Here, in an obvious notation, v denotes the right-handed neutrino superfield(s), ktH u v 
and k H \ LD are dimensionless coefficients, and we suppress flavor indices. The first 
term fl2.26ajl leads to Dirac neutrino masses when X attains its F-term VEV, (Fx) ~ 
m 3 / 2 Mp, while in the case of (12. 26b j) one has to observe that, due to the presence of the 
'non-perturbative' \i term, also H d attains an F term VEV, (F Hd ) ~ \i (H u ) ~ m 3 / 2 Uew- 
As qu u + qn d = mod M, both terms are allowed if q v + qn u + qi = mod M, which is 
precisely the condition that an effective holomorphic Y v term is allowed. Altogether we 
find, analogous to what we have discussed around (12. 2p . that effective neutrino Yukawa 
couplings 

Yv Mp" Mp (2 ' 27) 

will arise. For m 3 / 2 in the multi-TeV range this can lead to realistic Dirac neutrino 
masses. If we are to connect the suppression of Y v to the smallness of the \x term, it 
is natural to assume that the neutrino Yukawa coupling is forbidden by the same R 
symmetry that also forbids fi. As discussed above, L H u v has to have R charge 0. 
Moreover, there will also be holomorphic contributions to the Yukawa coupling. That 
is, even if both kiH u u and k H u- vanish, Dirac Yukawa couplings of the order m 3 / 2 /Mp 

d 

will get induced, where, as in our discussion of the fx term, m 3 / 2 represents the order 
parameter for R symmetry breaking. 



2.7 Discussion 

We have surveyed anomaly-free symmetries which forbid the \x term and are consistent 
with the Giudice-Masiero mechanism and SU(5). We find that these are discrete R 
symmetries Ti^ with M = 4m, m G N. The R charges of the H u are such that 
one expects a holomorphic contribution to the \x term of similar size. That is, the 
Giudice-Masiero mechanism strongly suggests the presence of additional holomorphic 
contributions to the effective \x term! 

Assuming further that the symmetries allow the up- and down-type Yukawa coup- 
lings and commute with flavor we find that they automatically forbid the troublesome 
dimension five proton decay operators and in many cases those of dimension four. In- 
terestingly, all these symmetries require a GS axion for anomaly cancellation. That is, 
these symmetries appear to be broken at the non-perturbative level. In other words, 
imposing compatibility with the Giudice-Masiero mechanism leads us to a situation in 
which a holomorphic \i term appears at the non-perturbative level, i.e. in a way the 
Giudice-Masiero term is unnecessary. 



3 Classification and models 

In this section, we explore anomaly-free discrete symmetries that solve some of the most 
severe problems of the MSSM. We will demand that the symmetry 

1. is flavor-universal and Abelian, i.e. a Zf f symmetry; 
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2. commutes with SU(5); 

3. forbids the /i term perturbatively; 

4. allows the usual Yukawa couplings; 

After revisiting in Section [37X1 the scan performed in [8], where Majorana neutrinos were 
considered, we turn to the Dirac case in Section 13.21 

3.1 Models with Majorana neutrinos 

In [SII2T], anomaly-free discrete R symmetries with qg = 1 were studied which satisfy 
the requirements [TH1] and in addition 

5. allow the Weinberg neutrino mass operator. 

It was found that there are only five phenomenologically attractive symmetries that 
commute with SU(5), one of which, a simple 7L^ symmetry, commutes also with SO(10). 
Further, the \x term, while perturbatively forbidden, appears at the non-perturbative 
level in four out of the five symmetries, and thus can explain its suppression. There is 
one symmetry which is anomaly-free without GS contribution; here anomaly freedom 
requires the number of generations to be a multiple of 3 [22] (for a similar connection 
between the number of generations and anomaly-free non-i? symmetries see [6|l30]). 

In the classification of [8], appears to be particularly attractive. Apart from the 
fact that it is the unique solution that commutes with SO(10), only lif provides a real 
solution to the \i problem. In this case, the discrete charges of H u and add up to 
mod M = 4 such that the \i parameter will be of the order of the gravitino mass, 
i.e. the order parameter of R breaking. This feature is not shared by the other four lifj 
symmetries, as also can be seen from our analysis in Section 12.41 In particular, it was 
argued that /i ~ m 3 / 2 for the case of lif . To substantiate these claims, an explicit string 
model with exact MSSM spectrum and the 7L^ symmetry was constructed in which the 
relation \x ~ (W) ~ m 3 / 2 is due to gauge invariance in extra dimensions |31j . 

Assuming in addition a Giudice-Masiero-like mechanism, one can see that TL^ is the 
unique solution also for general qg as follows. From the requirement that the Weinberg 
operator be allowed we infer that 



2q e mod M rx = q e — q Hu mod M/2 . 



(3.1) 



On the other hand, from the down-type Yukawa coupling it follows 



9io = - Ob ~ QH d + 2qe mod M 



f33J 



<?e + Qh u - qn d mod M/2 . 



(3.2) 



Demanding that the up-type Yukawa coupling be allowed leads to 



2q e - 2q 10 



mod M 



JB~2l 



2Qh u + 2q Hd 



mod M 



mod M , 



(3.3) 
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such that 5qn u = mod M. This means that qu u = mod M unless the order is a 
multiple of 5. In the latter case we can write the symmetry as x 22j^/ 5 where the 
factor is a non-i? symmetry. Hence we can focus on qu u = mod M, which implies, 
by fl2~3dl) . that q Hd = mod M. Then Equations flU]) and flX2]) imply 



qio = Q§ = qe mod M . (3.4) 

That is, the symmetry commutes with SO(10) in the matter sector. We already know 
from our discussion in Section 12.41 that the only meaningful R symmetry with this 
property is Tif. 

We also scanned the discrete Zf f symmetries up to order 200 with general qg without 
assuming a Giudice-Masiero-like mechanism. We obtain, apart from the symmetries 
of Tables 2.1 and 2.2 of [8], only a few new symmetries. However, as we show in the 
following in an example, these additional symmetries are redundant: consider a 
symmetry with (qio,q5,qH u ,qH d ,qe) = (1)17,8,52,5). This is equivalent to a Wif x Z 5 
symmetry with charge assignment ((1,3), (1, 1), (0,4), (0, 1), (1,0)). The Wi 5 is nothing 
but the non-trivial center of SU(5), i.e. it does not forbid any couplings (see the dis- 
cussion in [ISJE2]) an d the (non-trivial) factor is the one just discussed in the last 
paragraph. 



3.2 Models with Dirac neutrinos 

By modifying the above conditions, i.e. by demanding that the symmetry 

5. forbids the Weinberg neutrino mass operator perturbatively 
and 



6. is compatible with the Giudice-Masiero mechanism 

we obtain further interesting discrete R symmetries. Some sample symmetries are 
listed in Table [TJ Anomaly-free (non-i?) symmetries which allow for Dirac neu- 
trino Yukawa couplings have been discussed in [33]. The symmetries of Table [1] are 
inequivalent. One way of verifying this is to check whether or not two given charge 
assignments are equivalent by computing their Hilbert superpotential basis [34]. Only 
if the bases coincide, the assignments are equivalent. In the case of R symmetries, the 
Hilbert superpotential basis comprises homogeneous and inhomogeneous elements, or 
monomials. Every possible superpotential term contains precisely one inhomogeneous 
monomial and an arbitrary number of homogeneous monomials. In appendix O we list 
the Hilbert superpotential basis for examples with the !Zf 2 symmetries. 

3.2.1 Comments on the 7L^ symmetry 

One of simplest charge assignments appears to be the one of the 7L^ symmetry. Clearly 
the usual Yukawa couplings 10 10 H u and 10 5 Hd are allowed. Further, the Higgs 
bilinear H u Hd has R charge mod 8. If we assign the right-handed neutrino v R 
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(a) symmetries. (b) Residual symmetries. 
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Table 1: Classification of anomaly-free discrete R symmetries that forbid neutrino 
masses perturbatively. We restrict to orders < 36. (a) shows some sample symmetries. 
The equality between qo and p is due to Equation (I2.23p . The charge of the right-handed 
neutrino superfield v is determined by the requirement that q„ + qn u + qi = mod M 
(cf. the discussion below (12.26P ). In (b) we display the residual symmetries that remain 
after the ('hidden sector') superpotential acquires its VEV. 

charge 1, the Dirac neutrino Yukawa coupling will also be induced by R breaking. That 
is, we will have an effective superpotential which is schematically of the form 

W cS ~ 7713/2 H u H d + LH u u + QQQL . (3.5) 

IVlp IVl-p 

Here we suppress flavor indices. Once the superpotential of the hidden sector acquires a 
VEV, the lig is spontaneously broken down to a symmetry under which all matter 
fields have charge 1 and the Higgs fields have charge 2 (Tabled] (b)). Of course, this 
symmetry gets broken down to the usual matter (or 'i?') parity once the Higgs scalars 
attain their VEVs. 

The Hilbert superpotential basis [M] for this model (setting all quarks to zero) is 
given by the inhomogeneous monomials 

z/ 4 ; (LLE)u; (L H d E) ; (ill) 4 ; (L LE) 2 (L H u f ; (LH u f , (3.6) 



12 



while the homogeneous monomials are 

v 8 ; (LH u )u; (L H u f ; (L LE) 5 u ; (L LE) 4 (L H d E) ; 
H u H d ; (LLE) z/ 5 ; (L H d E) z/ 4 ; (L LE) 2 (L H d E) (L H u f ; 
(LLE) 8 ; (LH d E) 2 ; (LLE~) (LH d E~) u ; (L LE) 2 v 2 ; 
(LLE) 3 (LH u ) ; (Lff^) (L# u ) 4 ; [LLE) (LH u f . (3.7) 

Furthermore, there will be Kahler potential terms 

* D ^ {^ HuHd+ ^ LHu,y + ^ LQQQL ) +h - C - (3 ' 8) 

with X denoting the field that breaks supersymmetry, kn u H d , kLH u u and I^qqql being 
coefficients (and the flavor indices again are suppressed). The ku u H d term is nothing but 
the famous Giudice-Masiero term [7]. 

An important feature of this setting is that lepton number is violated at the quartic 
level, but bilinear lepton number violating terms are absent. That is, this model predicts 
the absence of neutrinoless double /3 decays. On the other hand, lepton number is not 
a good symmetry, which might have, for instance, important implications for the early 
universe. 

Let us also note that the coefficients in the above Kahler potential are not necessarily 
of order unity. In specific string constructions, these coefficients can in fact be as large as 
0(10—100) due to the presence of copious heavy states and/or combinatorical factors (cf. 
the discussion in [33]), enabling realistic predictions for neutrino masses in the sub-eV 
range. 

3.2.2 Comments on the symmetries 

Both Tif symmetries of Table [1] are problematic as they allow some R parity violating 
couplings. In particular, the first 7L^ allows for bi-linear R parity violation, i.e. the 5 H u 
coupling, while the second Wif admits the tri-linear R parity violating terms 10 5 5. In 
addition, both settings allow for a non-perturbative neutrino bilinear v v. That is, these 
symmetries can give us a non-perturbative Majorana neutrino mass term, which might 
be relevant for the construction of models realizing a TeV-scale see-saw scenario. Given 
our previous discussion, a straightforward possibility of rectifying this is to amend the 
settings by the residual symmetry from above (Table [2]). 

The Tif symmetries originally give us two inequivalent Hilbert superpotential bases, 
however, amending the settings by the above-mentioned TL^ symmetry leads to the same 
basis. Therefore, both 7L^ x symmetries give us the inhomogeneous monomials 

(LH d E) ; (LLE) u; (LLE) (LH u f ; (L LE)\L H u ) , (3.9) 
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(a) First Zf . 



(b) Second Zf . 



gio ffs gff d P gg 



9io gs Qh u qn d qe P qv 



1% 2 2 1 1 2 
Z 4 1 1 2 2 1 



Z£ 2 2 2 2 1 1 
Z 4 1 1 2 2 1 



Table 2 



Z4 extensions of the symmetries of Table [TJ 



whereas the homogeneous ones are given by 



v A ; H u H d ; z7 ; (L H u f ; (L H d E) (L H u ) (L LE) 3 

(LH d E) 2 ■ (LLE) (LH d E) v ; (L LE) 2 v 2 ; 
(LLE) 2 {LH u f ■ (LH d E) (L H u f (LLE) ; (L LE) A . 



(3.10) 



As before in our setting, bilinear lepton number violating terms are absent. In both 
cases this feature is due to the (anomaly-free non-R) symmetry, which commutes 
with SO(10) for the matter fields and is a consistent symmetry of the MSSM. Unlike the 
R symmetries, this symmetry does not forbid the \x term nor the dimension five proton 
decay operators. 

4 Summary 

The MSSM provides a very attractive scheme for physics beyond the standard model. 
However, in order to address its shortcomings, one, arguably, has to impose additional 
symmetries. Motivated by the structure of matter and the attractive picture of gauge 
unification, we have considered symmetries that commute with SU(5) in the matter sec- 
tor. From the requirement of anomaly freedom it follows that only discrete R symmetries 
can forbid the /x term. We also pointed out that anomaly matching for R symmetries 
in SU(5) symmetric models implies the existence of split multiplets below the GUT 
scale, with the simplest option being that a pair of Higgs doublets cancels the anomaly 
mismatch between the gauginos. Further demanding that a fi term of the order of the 
gravitino mass arises from supersymmetry breaking, i.e. either from the Kahler potential 
or from the non-trivial superpotential VEV in the 'hidden sector', we showed that the 
Higgs bilinear H u H d has to carry trivial R charge. We find that discrete R symmetries 
with these properties automatically forbid dimension-five proton decay operators once 
the usual Yukawa couplings are allowed. Even more, all symmetries appear anomalous 
such that a holomorphic \x term gets induced at the non-perturbative level. That is, de- 
manding compatibility with the Giudice-Masiero mechanism brings us to the situation 
in which a \x term of the desired magnitude appears even without the Giudice-Masiero 
term in the Kahler potential. 

We then discussed neutrino masses in the emerging MSSM models amended by dis- 
crete R symmetries. Restricting ourselves to flavor-universal Abelian, i.e. 52 jj^-, sym- 
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metries we find that, by demanding that the Weinberg operator LH U LH U be allowed, 
there exists only one possible symmetry, namely a 1^ symmetry Following a different 
approach, this Wif has also recently been shown to be the unique anomaly-free symmetry 
that commutes with SO(10) [2T]. The proof in [2TJ assumed that the charge of the su- 
perspace coordinate 9 can always be set 1, which we find to be too strong a requirement. 
However, we find that, if one is to allow for arbitrary 9 charges, this only leads to trivial 
extensions of Wif , such that the uniqueness of still prevails. 

If one requires instead the discrete symmetry to forbid the Weinberg operator, one 
can explain small Dirac neutrino masses. In particular, we successfully obtain a relation 
between the smallness of Dirac neutrino Yukawa couplings and the fi term which is 
based on anomaly-free discrete R symmetries with the above properties. Specifically, 
we find a class of anomaly-free discrete symmetries in which the appealing relations 
/i ~ (W)/Mp ~ m 3 / 2 and Y v ~ /x/Mp naturally emerge. 
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A Anomaly coefficients for symmetries with ar- 
bitrary qe 

The anomaly conditions for discrete R symmetries depend on qg. Consider a Zf f sym- 
metry, under which the superpotential transforms as 

W e 2 *^l M W (A.l) 

with q w = 2q e (such that jd 2 9 W is invariant). Superfields $ (/) = + y/29^ + 
99 transform as 

$(/) e WVM $(/) _ ( A _ 2 ) 
Correspondingly, the fermions transform as 

^(/) = e 2.i (g (/)- 9fl )/M^(/) (A3) 
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The anomaly coefficients hence read (cf. [H Appendix B], where the anomaly coefficients 
for the special case q$ — 1 are shown) 

£(rW) • - q e ) + q e £(adj G) , (A.4a) 

/ 

J2(Q if) ) 2 dim(r^) • (qW - q g ) , (A.4b) 

/ 

-21 q e + q e ^ dim(adj G) + ^ dim(r (/) ) • - . (A.4c) 

G f 

Here denote the charges of the superfields, the charges of the corresponding 
fermions are shifted by q$, q^u) = — qe- In Equation (IA.4al) . £{r^) denotes the 
Dynkin index of representation r^> normalized to £(N) = \ for the fundamental rep- 
resentation N of SU(iV) and £(adj G) = 02(G) represents the contribution from the 
gauginos, i.e. £(adj SU(iV)) = N. The first and second terms on the right-hand side of 
Equation (!A.4cP represent the contributions from the gravitino and gauginos. 



Aid 



M 



R 



U(l)-U(l)-S 



A p 

■^grav— gray— %yi 



B Green— Schwarz anomaly cancellation and anomaly 
universality 

In this Appendix, we discuss the discrete Green-Schwarz (GS) anomaly cancellation 
mechanism, following [8]. We start by reviewing the GS mechanism for a continuous 
U(l) symmetry in IB. It In IB. 21 we discuss the discrete version while lR3l is dedicated to 
the discussion of anomaly universality. 



B.l Anomaly cancellation for 'anomalous U(l)' symmetries 

We start by discussing the mixed anomaly coefficients G — G — U(l) anom for a simple 
gauge group G. There will be an axion a which couples to the field strength of G via 

-^axion 3 -F h F b . (B.l) 

o 

A possible prefactor can be absorbed in the normalization of a, which we do not specify 
here. Consider now the gauge transformation 

^(/) eio^QiSm^j) 5 ( B 2 ) 

where ip^ (1 < / < F) denotes the fermions of the theory and Qanom their charges. 
The crucial property of the axion a is that it shifts under (1B.2P as 

a — > a + 2^gs oc(x) . (B.3) 
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We can now fix the the Green-Schwarz coefficient 5qs from the requirement of in- 
variance of the full quantum theory. It follows from flB.3j) that, under a U(l) anom trans- 



formation with parameter a, the axionic Lagrangean shifts by 

A^ axion = -^-5 GS F b F b . (B.4) 
lb 

The Green-Schwarz term Sqs can now be inferred by demanding that the transformation 
of the axion a cancels the anomalous variation of the path integral measure [36"P37]. The 
latter can be absorbed in a change of the Lagrangean 

A.£f an omaly = ^^2 Fb ^ b ^G-G-U(l)anom " ( B ' 5 ) 

The coefficient A is the anomaly coefficient, given by 



r 



where the sum runs over all irreducible (fermionic) representations of G, 
denotes the Dynkin index of and Qanom is the U(l) anom charge. 

The axion shift allows us to cancel the G — G — U(l) anom anomaly by demanding 
AJz? anomaly + Aj2f axion = 0. This fixes the Green-Schwarz constant to be 

2tt 2 5 gs = A G _ G _ u(1)anom . (B.7) 
B.2 Discrete Green-Schwarz mechanism 

The Green-Schwarz mechanism also works if we replace U(l) anom by a discrete %u- 111 
this case the parameter a is no longer continuous but a = =jjjp with some integer n. Of 
course, there is no gauge field associated with the TLm- The discussion then goes as in 
the previous subsection. The discrete Green-Schwarz constant is now defined in such a 
way that under the 7Lm transformation of fermions 

tfS) e-^V^ (B.8) 
the axion shifts according to 

a a + ^A GS , (B.9) 

where A G s is fixed only modulo rj, 

7tMA gs = Ag—g~% m mod//. (B.10) 

The anomaly coefficients can be obtained from Equation (IB. 61) by replacing the U(l) anom 
charges Qanom by the %m charges q^ m \ 
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B.3 Multiple gauge groups and "anomaly universality" 

Let us now discuss the case of multiple gauge groups Gi. In heterotic string models very 
often the U(l) anom , Win or W^ anomaly coefficients fulfill certain universality relations, 

A Gi _ Gi _ H = p (B.ll) 

for all i and in this section H denotes either U(l) anom , Wn or Wfj. We will refer to (IB. lip 
as "anomaly universality" . In a recent paper it has been pointed out correctly that 
this may not necessarily be the case in general. That is, the anomaly universality (IB. lip 
is not a direct consequence of GS anomaly cancellation. 

In detail, multiple gauge groups Gi in general do allow us to introduce different 
couplings Q of the axion a to the various field strengths, 

^axion D 5>^?. (B.12) 

i 

The requirement that under an H transformation the contribution from the path integral 
measure gets cancelled by the discrete shift of the axion then implies that 

2vr 2 Q(5 G s = Aci-Gt-H (B.13) 

for all i. That means that the q coefficients can be chosen in such a way that the 
transformation of the path integral measure gets cancelled for each Gi gauge factor 
separately. In particular, one finds (in agreement with [38]) that in general the mixed 
Ad-d-H do not need to be universal. 

However, if this was the case in a given model, one would spoil the beautiful picture 
of MSSM gauge coupling unification. Let us spell out the argument in some more detail. 
In supersymmetry, the Lagrangean (IB. 12|) implies that there are couplings between the 
superfield S which contains the axion, S\g =0 = s + ia, and the supersymmetric field 
strengths associated to the gauge group factors Gi, i.e. 

^axion D J2 j^ 9 | S W * )W(l)a ■ (B.14) 

Once the real part of S acquires a VEV this will give rise to a non-universal change 
of the gauge couplings unless the q coefficients are all equal for the SM gauge group 
factors Gi = SU(3)c, SU(2) L and U(l)y. That is, anomaly universality is also required 
in order not to spoil the beautiful picture of MSSM gauge coupling unification. 

Furthermore, there might be model-dependent reasons why the Ag^g^h can be 
universal, for instance if all Gi come from a (for instance grand unified) simple gauge 
group, as we assume in the main body of this paper. Then the term 

=^axion ^ a -^GUT^GUT (B.15) 

is obviously gauge invariant. Hence, the anomalies need to be universal at least at the 
GUT level, as discussed around Equation ( I2.10p . 
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How could this universality possibly be broken? One may now worry about additional 
terms of the form 



where the operator $gut furnishes a non-trivial GUT representation (such as a 24-plet 
of SU(5)) and the parentheses denote a non-trivial contraction of the group indices H 
However, at the GUT level, i.e. for a trivial $gut VEV, such a term can not cancel the 
transformation of the path integral measure by a shift transformation of the axion a. In 
other words, it is not allowed by the symmetries of the action if we require that a shifts. 
Hence, these operators can not break anomaly universality. 

However, there is a second possibility. In higher-dimensional, e.g. in orbifold GUT 
type, models there can be localized terms which do not respect the GUT symmetry. 
That is, in settings where the GUT symmetry is broken locally in some regions of 
compact space such as orbifold fixed points, anomaly non-universality can arise. After 
integrating over compact space in order to derive the four-dimensional effective action 
one can indeed arrive at non-universal couplings c, of the axion to the three F F terms of 
the standard model. Still, as discussed before, if one is not to spoil the beautiful picture 
of MSSM gauge coupling unification, the A G ._ G ._ H coefficients need to be universal 
and these localized contributions have to be avoided. One possibility to avoid them 
is "non-local GUT breaking" in extra dimensions, which has been argued to yield the 
most appealing scenarios of precision gauge unification [^0l - H2] . In such scenarios, the 
localized dangerous GUT-breaking operators do not exist and hence the anomalies are 
universal @ 

Let us also comment on another statement in First, we would like to point 

out that the number of axions is not related to anomaly universality. Specifically, in the 
presence of multiple axions, which are available in heterotic compactifications [H2JII], one 
would have to define how they transform under a U(l) anom (or discrete) transformation. 
Since there is only one such transformation, this allows us to identify one unique linear 
combination of axions, called a as in our discussion above, which shifts while the other 
'would-be axions' stay inert. Therefore, the number of axions is not related to the 
question of anomaly (non-) universality. 

Furthermore, the authors of [38] argue that the anomalies cannot be universal both 
before and after doublet-triplet splitting. We disagree with this statement. First of 
all, 'before doublet-triplet splitting', i.e. before GUT breaking, there are more states 

3 The relative coefficients Cj of the axion coupling to the three F{ Fi terms of the standard model 
originating from <£>gut can be inferred from [39] . 

4 In Abelian orbifold models such operators can only stem from localized fluxes, which are Abelian 
(i.e. U(l)) fluxes. Hence, the A^-d-H coefficients coincide for all non-Abelian factors G{ from each 
Eg in such models. This is also in agreement with (38] where it is found that, in compactifications of the 
heterotic Es x Es string on blown-up orbifolds with Abelian fluxes, non-Abelian anomalies of each Eg 
factor are still universal. Since the relevant assumption in Section [2.21 needed to prove the uniqueness of 
7Lf (see also Section l2.2.6[) is that ^4su(2) L -su(2) L -z£ f an d ^.su(3)c-SU(3) c -Z R coincide, the uniqueness 
of Wif is also given in such constructions. 




(B.16) 
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around which contribute to the anomalies and anomaly universality follows from gauge 
invariance under the GUT group, see Equation (IB. 15}) . Moreover, in the absence of 
localized GUT breaking terms, if the anomaly coefficients are universal at the GUT 
level, where the contributions of extra states have to be taken into account, they should 
also be so in the MSSM. This is, again, nothing but 't Hooft anomaly matching (see 
Section |2.2.3[) . In fact, at the GUT level there is just one (unified) gauge group, such 
that universality is trivial. 

C The Hilbert superpotential bases for models with 
symmetries 

In Section 13.21 we discuss several TL^ symmetries that forbid neutrino masses perturba- 
tively and also present the Hilbert superpotential basis for a model with a Ti% symmetry 
and two Wif symmetries amended by an extra TL^ factor. In this appendix we provide 
further examples based on the TL^ 2 symmetries. As we have already stated above, ev- 
ery possible superpotential term ^ contains only one inhomogeneous monomial and an 
arbitrary combination of homogeneous monomials [M], i.e. 

^ = n with Vj g n , (c.i) 

where is an inhomogeneous and ^jom a homogeneous monomial. 

In Section 13.21 we list three examples which have a 1if 2 symmetry. As we will see in 
the following, the three sets of monomials differ. Hence, the three HJ^ 2 symmetries are 
inequivalent. The first symmetry has the charge assignment 

( ?io % Qh u q Hd qe P q, ) = ( 1 9 4 8 3 3 11 ) , (C.2) 
which leads to the inhomogeneous monomials 

(LH d E) ; (LH u f ; P 6 ; (LLE) P; 

(LLE) 6 ; (LLE) 4 (L H u f ; (L LE) 2 {LH u f , (C.3) 

whereas the homogeneous ones are given by 

(LLE) 12 ■ (LH U ) 12 ; H u H d ; (LH U ) P ; (LLE) (L H d E) P ; 
P 12 ; (LH d E) (LH u f ; (L H d E) 2 ; (LLE) 1 v ; 
(LLE) (LH u f ; (LLE) 6 (LH d ~E) ; (L H d E) P 6 ; 
(LLE) 2 P 2 ; (LH d E) (LLE) 4 (LH u f ; (LLE) 5 (L H u ) ; 
(LH d E) (LLE) 2 (LH u f ; (LLE) P 7 ; (LLE) 3 (LH u f . (C.4) 
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The second !Zf 2 symmetry has the charges 

( 9io % q Hu q Hd qe P qv ) = ( 2 6 2 10 3 3 4 ) , (C.5) 

which gives us for the inhomogeneous monomials 

(LH d E) ; (LLE) 3 ; (L LE) u; (L LE) (LH u f , (C.6) 

and for the homogeneous monomials 

(LLE) 6 ; (LLE) 4 9; (L H d E~) (L LE) 3 ; H u H d ; 
(LH U ) v ■ v 3 ; (LH u j 3 ; (LH d E) (LH u f (LLE) ; 
(LLE) 2 (LH u ) ; (LH d E) 2 ; (LLE) 2 v 2 ■ (L H d E) (L LE) v . (C.7) 

The last Zf 2 symmetry has 

( ?io % Qh u q Hd qe P q, ) = ( 4 10 2 3 3 2 ) , (C.8) 

as its charge assignment, with these we get the inhomogeneous monomials 

z/ 3 ; (LH d E) ; (LLE) v ; (L# u ) 3 ; (LLE) 2 (LH U ) , (C.9) 

and the following homogeneous ones 

{LH u f ■ (LH d E) (LH U ) 3 ; H u H d ; (L H d E) (LLE) 2 (LH U ) ; 

(LH U ) P; (LLE~) 3 ; (LLE) (LH u f ; (Ltf^) 2 ; 
(LH d E) v 3 ; (LH d E) (L LE) v ; (L LE) v A ; v & ; 

(LLE) 2 u 2 . (CIO) 
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